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1. Introduction 

What we observe visually as a single lightning flash actually has a 
complicated internal structure. The initial portion of a flash is called a 
stepped leader. The leader is forme'd when the cloud potential becomes 
sufficiently large to cause an electrical breakdown in the atmosphere. In a 
cloud-to- ground flash, the breakdown progresses in steps until it reaches the 
ground. When the leader reaches the ^ound a large current flows in the 
ionized channel caused by the breakdown and is called a return stroke, A 
sequence of several leaders and return strokes can follow. An entire 
sequence of associated leaders and return strokes is called a flash. A compre- 
hensive statistical model for the observed electromagnetic radiation generated 
by lightning must describe (1) the large scale flash structure, that is, the 
times between flashes and the amplitudes and durations of flashes and (2) the 
internal structure of a single flash including the portion generated by leaders 
and the portion generated by return strokes. 



The major effort during this reporting period has been devoted to modeling 
the large scale flash structure. Large scale flash data has been'measure'd’ 
manually from strip charts of storms of August 5, August 26, and Septerriber~12, 

1975. The data is being processed by a computer program called SASEV to 
estimate the large scale flash statistics.* The program, experimental'results, 
and conclusions for the large scale flash structure are described in Sections 2 

-and 3.- . z. s 

Some progress has been made in examining the internal flash's'ffucture" 

This consists mainly of developing the software required to process the NASA 
digital tape data on the University of Maryland’s UNIVAC 1108 computer and 
is described in Section 4. 

Section 5 lists plans for future research. 

2, Description of the Point Process Analysis Program SASEV 

P. A.W. Lewis, A.M. Katcher, and A. H.' IV’eis have’wfifteiTsrFORTKAN" ” " 
program called SASEIV for the statistical analysis of series ofevents [T]". The 
program implements the techniques described in the book. The Statistical 

Analysis of Series of Events , by D. R. Cox and P. A. W, Lewis -[--2]-: — We-have 

modified this program to run on the University of Maryland UNIVAC 1108 
computer in the batch or demand modes and have called the modified program 
SASEV. The principal modifications are the reduction of the length of the sequence 
of events that can be handled from 1999 to 1024 and the elimination of double 
precision arithmetic. Otherwise, the capabilities of SASEV and SASEIV are 
identical and are described in detail in Reference 1. The statistics computed 


2 



and tests performed by SASEV that have been found to be particularly useful 
in the analysis of lightning data are summarized in this section. 

2.1. Sample Moments 

Let . . . , be a sequence of N positive numbers. Normally, these 
'numbers would be the times between events in a point process. SASEV computes 
the following sample moments and normalized moments. 


(a) Sample mean 

, N 

^ J=1 ^n 


(b) Sample variance 

*2 1 y , " ,2 

n=l 


(c) Sample standard deviation 




-2 

(J 


(d) Coefficient of variation 


c = c/il 


(e) Third.central moment 




N 


(N-l)(N-2) n=l 


N 

2 (X -U)' 
n 


(f) Coefficient of skewness 

A A *3 

Y= U3/<l 
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(g) Fourth, cejitral moment 


^ N(N^ - 2N + 3) ^ _ ^^4 3(N-l)(2N-3) 

(N-l)(N-2)(N-3) n=l ^ (N-l)(N-2)(N-3) 


a 


4 


(h) Coeffitient of Kurtosis 


k = 



3 


Many physically ohserved point processes are accurately modeled as 

homogeneous Poisson processes. A reasonable first step in the analysis of a 

point process is to see if the Poisson model applies. A simple check is to 

compare the sample moments defined above with the theoretical statistical 

moments. The times between events in a homogenous Poisson process are 

independent, identically distributed random variables with an exponential 

-Xx 

probability density function, say f{x) - Xe ^(x) where u(x) is the unit step 
function. In this case the ideal moments are 
{a) M = E {x} = X"^ 


(b) a = E f(X - [if } = X‘^ 


(c) a = 


(d) c = a/^ = 1 


(e) H 3 = E {(X - u)^] = 2X'^ 

(f) Y= ^ 3 /^^ = 2 

(g) = E { (X - u)^ } = 9 X'^ 

(h) k = -3=6 
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2.2. The Sample Distribution Function, Log Survivor Function, and Histogram 


SASEV computes the sample distribution function for the sequence fx } as 


Fj^(x) = 


number of x • s ^ x 
n 

N 


This is an estimate of the actual distribution function F(x) = P f X ^ x} . 

The function R(x) = 1 - F(x) is called the survivor function. SASEV 
estimates R(x) as 

It computes and plots the log survivor function 


Gj^(x) = In Rj^(x) 


For the Poisson process, G(x) = In R(x) = -Xx. Thus deviations from linearity 
in the plot of Gj^(x) give insight into how the observed processes differs from 
an ideal Poisson process. Another closely related function is the hazard function 


h(x) 


d 

dx 


G{x) 


f(x) 

. 1 - F(x) 


The hazard function can be estimated from the plot of the log survivor function. 
For the ideal Poisson process, h{x) = X . 

SASEV makes a direct estimate of the probability density function for the 
random variables by computing and plotting a histogram for the observations 
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2.3. Goodness -of “Fit Tests for the Poisson. Process 

SASEV performs several goodness-of-fit tests to check if the homogenous 
Poisson process hypothesis is statistically reasonable. All of these tests are 
based on the fact that if a homogeneous Poisson process is observed over the 
interval (0, T), the observed normalized arrival times 

i 

y. = t./T = 2 X /T , i = 1,...,N 

i 1 1 — 1 » 

n-1 


are the order statistics of a random sample of size N from a population uniformly 
distributed over (0,1). Tests based on this fact are called uniform conditional 
tests. These tests are particularly convenient for the Poss son hypothesis because 
the rate parameter X need not be estimated and no grouping of data is required. 

Let the sample distribution function for the normalized arrival times be 



(y) = 


number of y. ^ y 
1 

N 


For the homogeneous Poisson process hypothesis, the theoretical distribution 
function would be 


jo f o r y ^ 0 
^(y) ^ > y for 0 ^ y ^ 1 

1 for y S 1 


SASEV computes the following four statistics. 

(a) The One-Sided Kolmogorov-Smirnov Statistics 

KS+ = Dn= V^^sop [F^(y) - y]= ^ max [-~ - yj 
O^y^l l^i^N 

~ ^ ^ sup [ y - Fj^(y)] = >/W max [y. - ] 

0^y<l l<i^N 
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(b) The Two-Sided Kolmogorov-Smirnov Statistic 


” y ! " ] 

0 1 


(c) The Anderson-Darling Statistic 


WN2 = = N 

N 


.1 [ yy) ’ yr 

J y(i-y) 


= -N - — S { (2i-l) In y. + [2(N-i) +1] In (1-y^) J 
1=1 


Each, of these statistics gives a measure of the deviation of ^j^(y) from the 
hypothesized distribution F(y) = y. The Anders on -Darling statistic emphasizes 
deviations near 0 and 1 as a result of the factor y(l-y) in the denominator. 

These statistics are used by selecting a threshold z and rejecting the 
hypothesis that the process is a homogeneous Poisson process, if the 
threshold is exceeded. The threshold is selected so that the test has a desired 
level a, that is, so that 


P i statistic > z 1 Hq i = a 

Asymptotic formulas for the levels of the Kolmogorov-Smirnov statistics 

are derived in Kendall and Stuart [3]. The asymptotic formula for the two-sided 

Kolmogorov-Smirnov statistic is 

00 

lim P[D >zj= 2 Z (-1)^“^ exp(-2r^z^) 

r=l 


The approximation is satisfactory for N ^ 80. Values of this expression for a 
range of z are given in Table 3 of Appendix 1. A table of asymptotic levels for 
the Anderson-Darling statistic is given by Lewis [4]. 
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The power of a test is defined to be the probability that the null hypothesis 

Hq is rejected given that an alternative hypothesis is true. Tests for the 

null Poisson hypothesis based on the uniform conditional property are not very 

powerful against a variety of alternatives. For example, suppose that the 

normalized arrival times {y^} are equally spaced in the interval (0,1). Then 

Fj^(y) remains close to F(y) = y. Durbin [5, 6] has suggested a modification. 

. of the •uniform, conditional tests that gives a large increase in power over the 

imiform conditional tests for a broad class of alternatives. The first step in 

Durbin’s modification is to order the sequence of times between events {x } 

n 

to obtain the observed order statistics 


0 ^ X ^ X ^ ^ X 

(1) “ (2) • • * " (N) 


The next step is to compute the sequence 


w 




(N + 2-1) 


^( 1 ) 


for 1=1, ...» N-1 


N 


■'N 


N 


-N 


where 


N 

t = S X 
N - n 

n=l 


It can be shown that under the null Poisson hypothesis, the sequence {w.} has 
the same distributional properties as the sequence . Thus, the Kolmogorov- 
Smirnov and Anderson-Darling tests described above can be applied to [w ] . 
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2.4* Serial correlation coefficients 


SASEV computes the serial correlation coefficients 


Z (x -p)(x -n) 

- _ _JJ_ i=:l " 

T ■ N-j N 

2 (X - \i ) 

i=l 


' A 

for j = 1, 2, . . , , min (N/2, 100) where y, is the sample mean. This sequence is 
an estimate of the autocovariance function 


p. = cov (X.,X )/var (X ) 
J 1 i+j 1 


Under the hypothesis that the x^'s are uncorrelated, that is, p = 0 for j 7 ^ 0, 

A 

the p ^s are approximately normally distributed with zero mean and variance 
J 

(N-j) . This approximation is reasonable for N ^ 100 if the skewness is 

moderate. SASEV plots (N-j)^p as a function of j. 

A renewal process is a point process in which the times between events 
are independent identically distributed random variables. The homogeneous 
Poisson process is a special type of renewal process. A set of independent 
random variables are always uncorrelated. However, a set of uncorrelated 
random variables may or may not be independent. Thus, the serial correlation 
coefficients can indicate whether or not a renewal process model is appropriate. 
2. 5. Spectrum of the Intervals 

SASEV estimates the spectral density of the, sequence of times between 


events as 

S{w) = ~ {1 + 2 


M 

2 

n=l 


A 


P 2- 


n n 


cos 
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where a is the Parzen window, that is, 
n 


"l.6{ 


+ 6 r 1 ” 1 f 


.uJ 

* * l M > 

a = ^ 

2 ( 1 -. 

Inl 

■) 

M 


n 



for j n 1 — M/2 

for — ^ ^ I n 1 ^ M 
elsewhere 


It also computes the discrete Fourier transform 



N 

2 

n=l 


^ 2TT(n-l)k/N 
n 


for k=0, . . . , N-1 


and the periodogram 



I \ \ ^/N 


for k=0, . . . , N“1 


It can be shown that if the times between events are independent identically- 

distributed random variables, then the values of C, for 0 < k < N/2 are 

k 

as-ymptotically independent identically distributed random variables with the 
exponential density. Consequently, the goodness -of-fit tests described in 
Section 2. 3 can be applied to the periodogram to check the renewal process 
hypothesis. SASEV appies the Kolmogorov-Smirnov and Anderson-Darling 
uniform conditional tests to the periodogram. 

2, 6. Spectrum of the Point Process 

SASEV estimates the spectrum of the point process in the following way 
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Let 


D = - t^)/(N-l) 


and 


N 

A(k) = 2 exp [j B (t - t )/D] 

i=2 ^ 


2 

I(k) = 2 1 A(k) 1 /(N-1) 


for k = 1, 2, . . , , P. B is normally chosen as 2 tt/ (N-1). P should be chosen 
to be larger than N. Usually all the important features of the spectrum will 
be shown if P = 2N. SASEV smooths I(k) over 5,10, and 20 points using a 
rectangular window and plots the results. 

Theoretically, the spectrum for a homogeneous Poisson process would be 

flat. 

3. Experimental Results and Conclusions on the Large-Scale Flash Structure 
Experimental data for storms of August 5, August 26, and September 12, 
1975 has been supplied by Goddard Space Flight Center in the form of strip charts. 
The time scale on the charts is highly compressed so that an individual lightning 
flash has a length in the order of one centimeter. Thus the envelopes but not 
the fine structure of the flashes can be observed on the charts. We are manually 
measuring the time between flashes, the peak amplitude of the flashes and the 
durations of the flashes from the strip charts and punching the data on cards 
in the format required by SASEV. A signal is judged to be a flash if its envelope 
exceeds 5 strip chart amplitude units on channel 1 and it also appears in several 
channels. The storm of August 5 has been completely measured and processed 
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by SASEV. The storm of September 12 has been completely measured and punched 
on cards but not yet processed. The storm of August 26 is approximately 95% 
measured. Some of the results of SASEV for the storm of August 5 are presented 
in Tables 1 and 2 in Appendix 1 and in Figures I to 60 in Appendix 2. They are 
discussed in this section. 

3.1. Time Between Flashes for the Storm of August 5, 1975 

•The results of SASEV for the time between flashes are summarized in 
Tables 1(a) and 2(a) and Figures 1 to 22. By the time between two flashes, we 
mean the time from the beginning of one flash to the beginning of the next flash. 

Table 1(a) shows the computed sample moments defined in Section 2. 1. 

The coefficients of variation are reasonably close to the ideal value of unity 
for a homogeneous Poisson process. The coefficients of skewness, for the most 
part, are close to the ideal value of 2. 0 for the Poisson process. The coefficients 
of kurtosis vary moderately from the ideal value of 6. 0 for the Poisson process 
but are still ”in the ball park. " Looking at the means, we see that the rate of 
flashes increases from tape 1 to tape 4 and decreases in tape 5. Thus, as 
expected, the flash rate is small at the beginning of the storm, increases in the 
middle, and decreases at the end. Clearly, the homogeneous Poisson process 
hypothesis is not valid on a long time scale since trends are present. 

Histograms of the times between flashes for tapes 1 through 5 are shown in 
Figures 1, 6, 11, 15, and 19, respectively. The corresponding log survivor 
functions are shown in Figures 2, 7, 12, 16, and 20. The histograms for tapes 
1, 3, and 4 appear as though they could be reasonably fit with an exponential. 
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The histogram for tape 2 is light at the low end while the histogram for tape 5 
seems heavy at the low and high ends. The deviations from an exponential 
density can be clearly observed in the log survivor functions as deviations from 
linearity. For the most part, the log survivor functions are reasonably linear. 

Table 2(a) displays the computed Kolmogorov-Smirnov and Anderson- 
Darling statistics defined in Section 2. 3. The column labelled KS shows the 
two-sided Kolmogorov-Smirnov statistic for the uniform conditional test. From 
Table 3 it can be seen that, using KS, the homogeneous Poisson process hypothesis 
is accepted at about the 80% level for tape 1, the 27% level for tape 2, the 71% 
level for tape 3, and the 39% level for tape 4. The hypotiicsis is clearly 
rejected for tape 5. The column labelled DN shows the two-sided Kolmogorov- 
Smirnov statistic when Durbin’s modification is used. Using DN, the homogeneous 
Poisson process hypothesis is accepted at about the 3. 9% level for tape 1, the 
0. 012% level for tape 4, and tie 0. 06% level for tape 5. It is clearly rejected 
for tapes 2 and 3. The uniform conditional test with Durbin’s modification is 
picking up definite deviations from the Poisson hypothesis. We conjecture that 
these deviations may be caused by trends and quantization of the measured data. 

The serial correlation coefficients for tapes 1 through 5 are shown in 
Figures 3, 8, 13, 17, and 21, respectively. The corresponding spectra are 
shown in Figures 4, 9, 14, 18, and 22. These plots indicate that, except in 
tape 5, the times between flashes are reasonably uncorrelated. There is definite 
correlation out to a lag of about 50 in Figure 21 for tape 5. This correlation 
causes the low-frequency peak in the spectrum in Figure 22. We conjecture that 
the correlation was caused by a trend. 
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The Kolmogorov-Smirnov and Anderson-Darling uniform, conditional tests 
were performed on the periodograms of the sequences of times between flashes 
as described in Section 2. 5. The results are presented in the last two columns 
of Table 2{a). The hypothesis that the times between flashes are independent is 
clearly accepted for tapes 1 through 4. It is accepted at about only the 0. 05% 
level by the Kolmogorov-Smirnov test and at about only the 0. 07% level by the 
Anderson-Darling test for tape 5. This is consistent with the conclusions from 
the plots of the serial correlation coefficients. 

The spectra of the point processes for tapes 1 and 2 are shown in Figures 5 
and 10 respectively. The plots are very wild indicating that the smoothing in 
SASEV IS insufficient. The points marked by x’s are the most smoothed. 
Smoothing the plots further by eye, we conclude that the spectra are consistent 
with the hypothesis of a flat spectrum for the homogeneous Poisson case. These 
spectra were not computed for the other tapes since the computation time is long 
and no unusual results were expected on the basis of the other computations. 

Taken as a whole, the results from SASEV indicate that the sequence of 
flashes on the large scale can be modeled as a homogeneous Poisson process over 
a time period that is moderate with respect to the duration cf the storm. 

3. 2, Peak Amplitudes of the Flashes for the Storm of August 5, 1975 

The peak amplitudes of the flashes were recorded in strip chart amplitude 
units and processed by the appropriate subroutines of SASEV. 

The various sample moments for tapes 1 through 5 are summarized in 
Table 1(b). Notice that the coefficients of variation, skewness, and kurtosis differ 
significantly from the theoretical values for an exponential distribution. 
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The Kolmogorov-Smirnov and Anderson-Darling statistics are presented 
in Table 2(b). Notice that the Kolmogorov-Smirnov statistic KS or Anderson- 
Darling statistic WN2 does not reject the exponential distribution hypothesis 
significantly but that with Durbin's modification the Kolmogorov-Smirnov 

statistic DN strongly rejects this hypothesis. These results are interesting in 

% 

that they show the lack of power of the ordinary uniform conditional tests against 
some alternative hypothesis. 

Histograms, log survivor functions, serial correlation coefficients, and 
spectra of the amplitudes are shown in Figures 23 to 41 for tapes 1 through 5. 

The histograms and, consequently, the log survivor functions differ significantly 
from tape to tape. Therefore, we propose no specific distribution function for 
the amplitudes at this time. However, it is clear from the plots that the 
exponential distribution is not appropriate. Except for tape 2, the serial 
correlation coefficients, spectra, and goodness -of -fit tests applied to the 
periodograms indicate that the amplitudes are uncorrelated. 

3. 3. Durations of the Flashes for the Storm of August 5, 1975 

The sample moments for the durations of the flashes in tapes 1 through 5 
are summarized in Table 1(c), 

Histograms, log survivor functions, serial correlation coefficients, and 
spectra of the durations are shown in Figures 42 to 60. The histograms do not 
seem to be fit by any well known densities. The log survivor functions have 
significant "deviations from linearity which rules out the exponential density. In 
particular, the histograms seem to be quite large for the small durations. The 
uniform conditional test statistics with Durbin's modification given in Table 2(c) 
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also strongly reject the exponential hypothesis. The corielation coefficients, 
spectra, and goodness -of -fit tests applied to the periodogram accept the 
hypothesis that the durations are uncorrelated. 

4. Progress in Investigating the Internal Flash Structure 

The major effort to date in investigating the internal flash structure has 
been the development of software to process the NASA digital lightning tapes on 
the University of Maryland's UNIVAC 1108 computer. 

The FORTRAN program CONVERT shown in Appendix 3 converts the 
16 bit two's complement format words on the NASA digital tapes to 36 bit one's 
complement format UNIVAC 1108 integers and then to 36 bit UNIVAC 1108 
floating-point words. The converted data can be plotted on the line printer using 
the plot routines in SASEV. 

The FORTRAN program PEAK shown in Appendix 3 detects peaks in the 
sampled lightning data. The following three conditions must be satisfied before 
a point IS declared a peak (1) an amplitude threshold must be exceeded, (2) the 
point must be greater than the adjacent two points, i. e. , it is a local maximum, 
and (3) the difference between the point and at least one of the adjacent points 
must exceed a threshold. PEAK was applied to a record of approximately 6, 500 
points and the detected peaks were compared with the plotted data. For this 
record, PEAK detected the peaks with 100% accuracy. 

The data on the NASA digital tape was selected to display the internal 
structure of a flash. The signals on the tape consist of a senes of similarly 
shaped pulses with varying amplitudes and spacing. It is believed that the pulses 
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are caused by the individual steps in the leaders. No indication of return strokes 
could be found in the file examined. A very preliminary analysis of the intervals 
between pulses has been made using PEAK and SASEV. For the small segment of 
data processed, the results of SASEV indicated that the homogeneous Poisson 
process model may not be appropriate for the sequence of pulses in the leader 
signal. Much more data must be analyzed before any definite conclusions are 
drawn. 

5. Future Research 

Research ideas to be pursued during the remainder of this grant and in 
the future, if funds are available, are outlined in this section. 

(1) We will complete measuring the large scale flash data from the strip 
charts for the storms of August 26 and September 9, 1975 and process the data 
with SASEV. 

(2) The large scale flash statistics for new storm data supplied by NASA 
will be processed with SASEV. 

(3) The instantaneous flash rate, that is, the flashes per minute, will be 
estimated for the storms on which we have data and plotted as a function of 
time. A program to pe'rform this task is nearly completed. In addition to 
providing a history of the storm intensity for NASA, the results will indicate to 
us if the tests performed by SASEV were influenced by trends. 

(4) A significant effort will be devoted to investigating the internal flash 
structure. Statistical models for the times between pulses and amplitudes of 
the pulses in the leader process will be developed using the results of SASEV. 

The sequences of return strokes will be investigated similarly. 
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(6) If time and computer funds permit, we will perform a spectral 
analysis on the digitized storm data supplied by NASA. Spectral analysis can 
often uncover interesting signal characteristics. 

(7) We will investigate modeling the entire lightning signal by a branching 
renewal process. A branching renewal process is a point process in which each 
event in a primary sequence of events spawns a secondary sequence of events. 
The primary events would be the beginnings of flashes and the secondary events 
would be the pulses within a flash. Methods for estimating process parameters 
will be investigated. 

(8) The leader process contains a sequence of similarly shape pulses 
with random amplitudes and spacing. We believe that the shape of these pulses 
is the system impulse response. The pulses often overlap and are corrupted 
by noise. We would like to investigate the problem of optimally estimating the 
impulse response of a system driven by a continuous -time point process when 
the observed data is the system output corrupted by additive noise and sampled 
at a uniform rate. To our knowledge, this problem has not been previously 
solved. 
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APPENDIX 1 
TABLES 


20 



■rAJBJbiU 1 


Storm of August 5, 1975 
Some Statistics Generated by SASEV 


TAPE 

(a) Time Between Flashes (Seconds) 

COEFF. OF 

N MEAN STANDARD DEVIATION VARIATION 

|i a a/fi 

COEFF. OF 
SKEWNESS 

COEFF. OF 
KURTOSIS 
3 

1 

97 

3.4876 

3. 0484 

0.8741 

2. 0935 

8.9313 

2 

410 

2.9318 

2. 3470 

0.8005 

1. 3940- 

4.9723 

3 

506 

1. 7802 

1. 6425 

0.9226 

1. 9734 

7.9215 

4 

240 

1.4079 

1. 3991 

0.9937 

1.9350 

6.8081 

5 

228 

6. 1395 

6.8915 

1. 1225 

1.6741 

5. 6559 

(b) Peak Amplitude of Flashes (In Strip Chart Amplitude Units) 

1 

97 

14. 9897 

5. 0796 

0. 3389 

0. 6523 

2. 3001 

2 

410 

11. 7829 

5. 6624 

0.4806 

0.7579 

2.5375 

3 

506 

10. 3945 

6. 5425 

0.6294 

1. 3719 

3. 5116 

4 

240 

8. 8500 

5. 1823 

0.5856 

1. 9760 

5. 9478 

5 

228 

10. 2939 

6. 1319 

0. 5957 

1.0990 

, 2.9002 

(c) Duration of Flashes (Seconds) 

1 

1 

97 

0.4299 

0. 2732 

0. 6354 

0. 8323 

*2.8771 

2 

410 

0. 3624 

0. 2400 

0.6622 

0. 9291 

3. 3349 

5 

506 

0. 2544 

0. 2288 

0.8990 

1.9212 

7. 1981 

4 

240 

0.2033 

0. 2263 

1. 1128 

3. 3910 

16.9794 

5 

228 

0. 2272 

0.2001 

0.8809 

2. 0493 

7. 0931 
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storm of August 5, 1975 

Kolmogorov-Smirnov and Anders on -Darling Goodness -of-Fit Tests 
for Exponential Density Based on Uniform Conditional Property 


(a) Time Between Flashes WITH DURBIN'S APPLIED TO 

TAPE N KS WN2 TRANSFORMATION ’ PERIODOGRAM 

DN WN2 KS WN2 


1 

97 

0.6321029 

0. 383976 

1.401567 

3. 550858 

0. 8107939 

0.9660473 

2 

410 

1. 045668 

2. 098625 

2.639786 

16. 57763 

0. 6290504 

0. 4334488 

3 

506 

0. 6995085 

0. 7348518 

3.317009 

12.56082 

0. 3746214 

0. 1418037 

4 

240 

0. 8894655 

1. 157148 

2. 196101 

4.452011 

0. 5797793 

0. 2736702 

5 

228 

2. 877865 

26. 99329 

2.004672 

7. 86273 

2. 034177 

6. 590207 


(b) Peak Amplitude of Flashes 


1 

97 

0. 2439943 

0.0568504 

5.440877 

60. 25202 



2 

410 






1. 558849 

4.433937 

3 

506 

0. 7116409 

0. 7573547 

10.82038 


0. 8282011 

0.8978100 

4 

240 

0. 5034376 

0. 3427868 

8.734251 


0. 5788312 

0. 2426233 

5 

228 

0. 5596259 

0.4229412 

7.318207 


0. 5745411 

0. 3834915 


(c) Duration of Flashes 


1 

97 

0. 7276520 

0. 7565508 

2.622042 




2 

410 

0. 8409247 

1.455425 

5.579909 


! 0.7827463 

0.9495754 

3 

506 

1. 632041 

4.418961 

8.840851 

56. 30646 

0. 9280081 

1. 249138 

4 

240 

1. 015736 

0. 8823013 

7.603098 

48.96658 

0. 7042681 

0.7879419 

5 

228 

0.9044824 

0. 9611149 

6. 631598 

32.84324 

0. 6494161 

0. 5272541 
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TABLE 3 


Asymptotic Significance Levels for the 
Kolmogorov- Smirnov Statistic 



z 

P( I>N>z ) 




'"'"7-rooooooo+oi"''" " " 

1 

y. 

' . 15 

. 10000000+01 



-20 

. 10000000+01 

» 


.2b 

.99999999+00 



• 30 

.99999070+00 


t 

.35 

.99969714+00 



;40 

.99719233+00 

- 


./j5 

, - .98741063+00 

> 

* 

. 50 

.96394525+00 



. 55 

.92281682+00 

- \ 

- 

' . 60 

.86428282+00 

> 


.65 

.79201308+00'- 

» 


. 70 

.71 123525+00 



. 75 ' 

> .62716711+00 

* 


.80 

.1344.1 /.2/-8 + 00 



.85 ' 

. 46531929+00 

\ 

\ 

.90 

.39273078+00 


t 

.95 

. 32748555+00 



1 .00 

' ,26999973+00 

\ 

- 

1 .05 

.22020562+00 

ji 

K 


1.10 , 

. 17771325+00 

\ 

1 


1 . 10 

. 14195991 + 00 


' 

1.20 ‘ 

. 1 12249 71 + 00 

i 

t 

1 .25 

.87866448-01 

\ 


1 . 30 

' .68092250-01 

X 

* 

1 .35 

".5224 191 1-01 



1 .40 

.39681899-01 

* 

y 

» \ 

1 .45 

.29841489-01 

• 

V 

1.50 

.22217975-0! 

\ 

r 

1 

i. 

1 -.55 

. 1637 74 12-'01 

\ 

1 .60 

. 1 1952054-01 

\ 

X 

i .65 

.86356858-02 

\ 


1,. ?0 

.61774354-02 

\ 


r. 75 

.43749856-02 

\ 

1 

\ 

1 .80 

.30676239-02 



1 .85 

•21295343-02 

4 

\ 

1 .90 

. 1 4636061-02 

1 

1 

1 

i 

1 .95 

.9959 1 1 77-03 

1 

( 

» 


2.00 

. 67092590-03 

( 


2.05 

.44749213-03 ' 

i 

t 

i 


2.10 

.29549713-03 

4 

r 


2.15 

.19318713-03 

1 

► 

2.20 

. 12504322-03 

t 

1 


2.25 

.80130745-04 


OEIGINAL PAGE IS 

2.30 

.50838800-04 


OF POOR QUALiry 

2.35 

.31933641-04 



2.40 

. 19859059-04 



2.45 

. 12227169-04 



2.50 

.74533286-05 

( 

i 

. 2.55 

,■44981264-05 
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FIGURES 
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HISTOGRAM OF TIME 0E TWEEN EVENTS 
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FIG. 3 SERIAL CORRELATION OP INTERVALS (TAPE 1) 
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PIG, 5 SPECTRUM OP THE COUNTING PROCESS (TAPE 1) 
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PIG. 6 HISTOGRAM OF TIME BETWEEN EVENTS (TAPE 2 ) 
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33 



FIG. 8 SERIAL CORRELATION OP INTERVALS (TAPE 2) 
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FIG, 10 SPECTRUM OP THE COUNTING PROCESS (TAPE 2) 
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FIG. Ill SPECTRUM OF INTERVALS (TAPE 3) 
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PIG. 29 SPECTRUM OF AMPLITUDES (TAPE 2) 
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LOG SURVIVOR FUNCTlOK 


AMPLITUDE OF FLASH 
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FIG. 31 LOG SURVIVOR FUNCTION (TAPE 3 ) 
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FIG. 32 SERIAL CORRELATION OF AMPLITUDES (TAPE 3) 


PLOT OF SERIAL CORRELATION OF AMPLITUDES 
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FIG. 36 SERIAL CORRELATION OF AMPLITUDES (TAPE k) 
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FIG.ilO SERIAL CORRELATION OP AMPLITUDES (TAPE 5) 
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fig.41 spectrum op amplitudes (tape 5) 
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FIG. 48 SPECTRUM OF DURATIONS (TAPE 2) 
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FIG. 51 SERIAL CORRELATION OF DURATIONS (TAPE 3) 
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FIG. 59 SERIAL CORRELATION OF DURATIONS (TAPE 5) 
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APPENDIX 3 
PROGRAMS 
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SASEV 


A Program for the Statistical Analysis 
of Series of Events 

SASEV IS a large EORTEAN piogram (approximately 3, 000 cards) for 
the statistical analysis of point processes. The program is documented in 
detail in Reference 1. A listing of SASEV will be provided on request. 
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CONVERT 


A Program for Converting NASA Dxgital Lightning 
Tapes to Oniversxty of Maryland UNIVAC 1108 Format 


Clf'EiNSICN r'UA i C 1 ^900)^ lDUi» < 1 & > ^ Fi\ C 60C0 ) 

ECUl valence: {DAi£i,FDftl) 

DAi* IDijy/l&V ’/ 

INitCEF PA^1 lY^D£NS^PLCvr»Sl^' lUb^l UF ( b /OO ) / DA J A < 12900) 
DAjA NFiLL/G (^77 77 7600000/ 

S1A1US=0 
PANl ] r-= 1 
DEi\S= 1 600 
kCCv-P=0 
i\..ChDS=5689 
N$A,%p= 12800 

CALL aREAC( 10>e>.'F>MCRDS>S'lAiUS>PA,sJ i Y # DENS> FCGv{0 
L= 1 

DC 1 1=1,MCaDS>/< 

DA 1A(L) =FLDC0> 1 6» E‘JF( 1 ) ) 

DA 1 A(L+ 1 )=FLiJ( 1 o> 16^PUF( i ) ) 
i\ 1 1 =FLC( 32/ Pur < i ) ) 
tvlR=FLD(0/ 12/FuFCI+l )) 

DAlACL+2) = i\ i i .LS. 1?. V.NT2 
DAlA(Lt3) = FLD<: 12/ 16/PUF( 1+ 1)) - 

M 1 =FLD(28/&/ BUF( 1+ } ) ) 

M 2=FLD(0/ &/ 8UF< 1 + 2) ) 

0AiA(L + A)=Nil.LG.&W.i\12 
DA 1 A(L+b> = FLD(6/ I 6/ PUF C I + 2)') 

Ni 1=FL0(2A, l2/F'JF(i + 2) ) 
t\12=FLD(0 / A / 5UF< 1 + 3) ) 

DA 1 A(L+ 6) =i\ 1 1 .L V J 2 
DA 1 ACL+ 7)=FLD(^7 1 6/ EU=^'( 1 + 3) ) 

DA i AfL+&)=FLD(20/ I 6/BL'F( 1 + 3) ) 

L=L + y 

DC2 I = l/N3A,/p 

IF(FLDC20/ 1 / DA i AU ) ) .EC. 1 ) DA 1 A ( 1 ) = ( DA i A ( 1 ) ~ 1 ) . v . Nr 1 LL 
DC 3 i=l/.\SAi^f- 
FDAi ( 1 )=DA 1 AC I ) 

EI= *E 1 • 

GC-='6G‘ 

1 1 ab=LN'= 1 

UkI 1E(2/A) CC/ LN-/ I i AP, 1 DU^ 

FCRNAI(P;c,A 2/ 212/ IbA/) 
tknE(2/A) E f/Li\/ 11 AP, IDON 
END FILE 2 
REi IND 2 
N = 600 
M = I 

DC b I=l/N 
FN (l)=i-l+Nl 

call FLlCI/|/l/l/0/0/2/l,l/0,10/t\/10/10/0/0» 

1 l./0/0/0/0/2/f’DAi(i\))/Fi\/i\/0/ 

2 O/O/O/O/O/O/O/O/O/G/O/O/O/O/O/O/G/O/O/O/O/O) 




IS 





PEAK 


A Program for Detecting the Peaks in Sampled 
Lightning Data 


DIKEKSIC^ XC665C) 

DIKEl\SIC!\' A^H( 600 ) 
htuliND- 3 
hCADC3) X 

iA=lCO. _ . ' 

1D=60. 

N = 0 

DC b l=2^66-^'9 
I"~CaCJ).LE.1A) CC I£ 5 

IFf X( 1 - ! ) .GE .XU ) .C ) .LI .A( 1+ 1) ) LC iC 5 

IJI=A<1)-X{1-1) 

U2=a< 1 )-xa-n > 

Ir CUl «L r. . 10) CC 1 C 5 

N = i \+1 
ll=l+6lb0 
A<K) = 1J 

Ai?'iP<rX)=X < 1 ) ■' 

5 CCMIMJE 

, V, 

X FCKKAicah \’=,no) 

tFllE<^>2) (X(i)^i=l,N) 

EKD i L£ A 

2 FC^^ A i C 12F6.0) 

\ n I 1 E ( 7 ^ 3 ) { Aiv r ( J ) ^ 1 = 1 ^ i\ ) 

£l\D PILE 7 

3 FCRt^A 1 clSFb.O) 

, , t EXD . 


OUlGINAj;, Pagp ra 
0 ® iWE 
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